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ABSTRACT 

We present results from Chandra observations of SDSS J1004+4112, a strongly lensed quasar system 
with a maximum image separation of 15". All four bright images of the quasar, as well as resolved 
X-ray emission originating from the lensing cluster, are clearly detected. The emission from the 
lensing cluster extends out to approximately We measure the bolometric X-ray luminosity and 
temperature of the lensing cluster to be 4.7 x 10 44 ergs -1 and 6.4 keV, consistent with the luminosity- 
temperature relation for distant clusters. The mass estimated from the X-ray observation shows 
excellent agreement with the mass derived from gravitational lensing. The X-ray flux ratios of the 
quasar images differ markedly from the optical flux ratios, and the combined X-ray spectrum of the 
images possesses an unusually strong Fc Ka emission line, both of which are indicative of microlcnsing. 
Subject headings: galaxies: clusters: general — gravitational lensing — quasars: individual 
(SDSS J100434.91+411242.8) — X-rays: galaxies 



1. INTRODUCTION 

The quadruply lensed quasar SDSS J1004+4112, first 
identified as a gravitational lens candidate from opti- 
cal imaging and s pectroscopy from the Sloan Digital 
Sky Survey (SDSS: lYork et al.ll2000ft has an exception- 
ally large image s eparation of ~ 15" (jlnada et al.ll2003t 
lOeuri et al.ll2004|) . The z = 1.734 quasar is multiply im- 
aged by a cluster of galaxies (rather than a galaxy) at 
z = 0.68. This is the only known quasar strongly lensed 
by a central part of a massive cluster. The quasar ap- 
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pears to be radio-quiet since it is n ot detected in radi o 
surveys such as the FIRST survey ( Bec ker et aTH l995). 
Subsequent observatio ns revealed additiona l intriguing 
aspects of the system. iRichards et alJ l)2004l) discovered 
time variability in the blue wings of several optical broad 
emission lines, which can be plausiblly interpreted in 
terms of microlesning. H ubble Space Telesc ope images 
detected the fifth image l)Inada et al.| |20 05[) and mul - 
tiply imaged background galaxies ijSharon et al.ll2005[) . 
SDSS J1004+4112 appears to be an ideal laboratory for 
exploration of the structure of a large cluster of galaxies 
and a quasar. 

The uniqueness of SDSS J1004+4112 affords advan- 
tages that strongly argue for multiwavelength observa- 
tions, of which X-ray observations constitute an essential 
part. Measurement of the X-ray properties of the lens- 
ing cluster offers an opportunity to improve lens models 
of SDSS J1004+4112. Moreover, the lensing cluster is 
the first example of a strong-lens selected cluster; it is 
interesting to determine if the cluster follows the em- 
pirical scaling relations between luminosities, tempera- 
tures, and masses. Second, we can investigate the mag- 
nification ratios of the quasar images at X-ray wave- 
lengths, which may aid the interpretation of the mi- 
crolensi ng event that has bee n observed at optical wave- 
lengths ijRichards et, al.1 2004). Indeed, X-ray flux ratios 
of lensed quasars are frequently different from optical 
flux ratios; this phenomenon has been attributed to m i- 
crolensing (jChartas et al.ll2004t iBlackburne et alJl2006|) . 

In this paper, we report our analysis of an observation 
of SDSS J1004+4112 with the Chandra X-Ray Obser- 
vatory. A key feature of Chandra images is the high 
angular resolution (on the order of 1"); this resolution 
is required to allow reliable separation of the X-ray flux 
of the lensing cluster from that of the quadruply im- 
aged quasar. In this paper we will adopt a cosmological 
model with the matter density Am = 0.27, the cosmo- 
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logical constant = 0. 73, and the Hubble constant 
H = 70 kms^Mpc" 1 <|Spergel et all 12003). At the 
redshift of the cluster (z — 0.68), 1" corresponds to 
7.16 kpc. Unless otherwise specified, quoted errors in- 
dicate the 90% confidence range. 

2. OBSERVATIONS 

SDSS J1004+4112 was observed for 80 ks with the 
Chandra Advanced CCD Imaging Spectrometer (ACIS; 
iGarmire et"al] l2003j) on 2005 January 1 and 2. The 
data were obtained with the ACIS S3 CCD operating 
in VFAINT mode. This CCD has a 1024 x 1024 pixel 
format with an image scale of / /492 pixel -1 . The target 
was offset from the nominal aim point with a Y-offset of 
— 1', however, this has little effect on the spatial resolu- 
tion. The CCD temperature during the observations was 
-120°C. 

The data were processed using the standard software 
packages CIA0 3.2.1 and CALDB 3.0.3. In the analysis 
of diffuse emission from the cluster (33J), periods of high 
background rates (defined by > 20% higher than the qui- 
escent rates in the 2.5-7 keV band) were removed; the 
net exposure time is 62042 s. On the other hand, the 
data without lightcurve-filtering (79987 s) are used in the 
analysis of quasar components (2J since the background 
counting rate is negligibly small for point sources. 

The full band (0.5-7 keV) ACIS image is shown in 
Figure □ The four images of SDSS J1004+4112 are 
clearly detected and resolved. We found astrometry off- 
sets for the Chandra data, (ARA, ADec) = (0'/4 ± O'.'l, 
— O'.'l ± O'.'l), from a comparison of three sources that lie 
within 65" from the aimpoint with those identified by the 
SDSS survey. After correcting for the mean offsets, the 
measured coordinates of the four X-ray quas ar images are 
consi stent with those obtained with Subaru ijOeuri et alJ 
2004) within < 0'.'2. In addition, extended emission, 
roughly centered on the quasar images and extending out 
to r ~ 1.'5, is seen in the X-ray data. This component 
should arise from hot gas in the lensing cluster. The rel- 
ative offsets of the coordinates of image A and the peak 
of the diffuse emission are (ARA, ADec)=(-7'/4 ± l'.'O, 
4'.'7± l'.'O); the location of the peak of the extended emis- 
sion is consistent with the position of the brightest cluster 
galaxy Gl. 

In the ACIS-S3 field, 39 point-like sources including the 
four quasar images were detected with the wavdetect al- 
gorithm with a significance threshold parameter of 10~ 6 . 
The bright source at (10:04:34.262, +41:12:20.33) with 
an X-ray flux of - 1.0 x 10 -13 ergs^^m" 2 (0.5-7 keV) 
has been observed serendipitously as part of the follow- 
up program to spectroscopically monitor components A 
and B at APO (Richards et al. 2004); a low S/N spec- 
trum indicates a redshift of z ~ 1.26, which is marginally 
consistent with the second most likely photometric red- 
shift (z — 1.325 with 26% confidence). The power- 
law slope, a ox , connecting the rest frame 2500A and 
2 keV flux densities, is ~ —1-2, which is roughly con- 
sistent with th e q n v -UV luminos ity re lation of AGNs 
l)Strateva et alJl2003 Isteffen et alJl2006j) . The source at 
(10:04:33.653, +41:13:07.24) has an X-ray counting rate 
of (4.4 ± 0.7) x 10~ 4 counts s , and it is not a counter 
lens image of the southern bright source as the separation 
angle is too large. 



3. LENSING CLUSTER 
3.1. Spectral Analysis 

We derive the spectrum of the cluster component by 
extracting the data from the circular region within a ra- 
dius of 1' that is centered on the position of the galaxy 
Gl. The fluxes from the four quasar images and several 
additional point sources were removed from the spectral 
integration area by excluding all regions within a radius 
10 times the size of the point spread function (PSF) at 
the source positions. Here the size of PSF is defined as 
the 40% encircled-energy radius at 1.5 keV at the source 
position. The background was estimated from the data 
in a surrounding annulus (2' < r < 2.'5). The source 
and background counts within r < 1' are 1629 ± 54 and 
906 ± 20 respectively. The 0.5-7 keV extracted spectrum 
is fitted with the MEKAL thin-thermal plasma model 
llMewe et al J 119851 119861 iKaastral Il992l iLiedahl et alJ 
|1995|) utilizing the XSPEC version 11. In this exer- 
cise, the Galactic hydr ogen column density is fixed to 
A H = 1.13 x 10 20 cm- 2 (|Dickev fc Loekmanl H9901. The 
X-ray temperature is constrained to be kT — 6.4l 2 ' 3 keV 
and the metal abundance to be Z — 0.21Z Q (the 90% up- 
per limit is 0.62Z Q ). The Galactic absorption-corrected, 
0.5-7 keV flux is 1.62 x 10 -13 ergs^cm- 2 (r < 1'). The 
bolometric X-ray luminosity within r§oo is estimated as 
Lx = 4.7 x 10 44 ergs -1 , where r 50 o is defined as the ra- 
dius within which the average matter density is equal to 
A c = 500 times the critical density of the Universe at 
the cluster redshift. The /3-model analysis in S 13. 21 yields 
?"500 = 0-79 Mpc. This luminosity is lower than the mean 
value expected from the luminosity-temperatu re relation 
of dis tant clusters, L x = Ot^xlO 45 ergs -1 (|Ota et alJ 
2006), but our estimated value is within the observed 
scatter of the data. The reduced chi-square of the best- 
fit model is x 2 / d of = 31.4/42. 

To investigate the radial temperature profile, we fur- 
ther analyzed spectra integrated from the inner r < 
0.'23(- 100 kpc) region and the outer 0.'23 < r < 1' re- 
gion of the lensing cluster. The fitted temperatures are 
kT = 6.3tl 5 7 keV and 5.9±?J keV, respectively. These 
results suggest that there is not a strong temperature 
drop towards the cluster center, indicating radiative cool- 
ing is not important in this cluster. 

3.2. X-ray Surface Brightness Profile 

The one-dimensional radial surface-brightness profile 
of the extended X-ray emission was constructed by 
adopting a center at Gl and azimuthally averaging the 
0.5-5 keV image, which was corrected for the tele- 
scope's vignetting and the detector responses and re- 
binned by a factor of two (i.e., 1 processed image pixel 
is 0'.'98). As shown in Figure [3 the radial profile was 
fit with the following two models: (1) the conventional 
/3-model S(r) = S„[l + (r/r c ) 2 ]- 3 ^ +1 / 2 , and (2) the 
pr ofile derived from the universal mass profile proposed 
by iNjiy^n^^^alJ 1^ 9971) plus isothermality of the clus- 
ter l|Suto et alJl!998l hereafter NFW-SSM). The back- 
ground was set to a constant in the fitting. We find that 
both models can fit the observed radial profile reason- 
ably well: The reduced chi-square is x 2 /dof = 213.2/196 
and 205.5/196 for the /3-model and NFW-SSM model, re- 
spectively. The fitted values of the /3-model parameters 
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are = 0.59±g;g| and r c = 9.Ttl° arcsec(= 69±J| kpc), 
while the NFW-SSM model results in B = 9.1±J:| and 
r s = 39tg 2 arcsec(= 278±f| kpc). The prohle of the 
intra-cluster medium (ICM) has a compact core, typical 
of observed surface brightness profiles of relaxed clusters. 

The fit shows excess flux in the measurements com- 
pared to the models within r < 15 kpc (Figure |2J- 
This is o ften seen in central regions of r elaxed clus- 
ters (e.g., iMohr. Mathiesen fc Evrardl 11999ft . The cen- 
tral emission is significantly extended compared with the 
PSF. The 3-cr upper limit on a point-source luminosity 
is ~ 4.4 x 10 42 ergs -1 (0.5-5 keV) assuming a power-law 
spectrum with index r = 2 and a source redshift of z = 
0.68. The excess luminosity within r < 2" is estimated as 
~ 7 x 10 42 ergs -1 (0.3-8 keV), which is higher than for 
a typical elliptical galaxy and rather c omparable to that 
of a cD galaxy in nearby clu sters fe.g.. lMatsushita et al.l 
120001 lO'Sullivan et aljfeoOl^ . Note that we subtracted 
the best-fit /3-model component from the total emission 
assuming the MEKAL model with kT — 6.4 keV and 
Z = 0.21Zq for the cluster spectrum. The B-band lu- 
minosity of the galaxy Gl is Lb ~ 3 x 1 n .Lm from the 
g magnitude of 22.11 ijlnada et alJl2003[) and the color- 
trans formation law for elliptical galaxies l)Fukugita et alJ 
1995). The above values may f ollow the Lx — Lb re- 
lation for elliptical galaxies (e.g.. lO'Sullivan et al]l200ij) 
and the excess flux may be attributed to hot ISM emis- 
sion from the central galaxy. However, the X-ray tem- 
perature inferred from the hardness ratio is 2.3^f' keV 
assuming a, Z = IZq MEKAL model. This is marginally 
high for a typical elliptical galaxy. Thus we might be 
looking at a superposition of the galaxy component and 
central emission from the cluster. Given the present pho- 
ton statistics, however, it is not possible to constrain 
further the origin of the emission. Note that the to- 
tal X-ray luminosity of unresolved low-mass X-ray bina- 
ries is expected to be of the order ~ 10 41 ergs -1 (e.g., 
iKim fc Fabbiand 12004) and contributes only ~ 1% to 
the excess l uminosity. X-ray emission from the central 
fifth image ijlnada et al.ll2005[) is unlikely: the contribu- 
tion is estimated to be only ~ 3 — 9 counts (i.e., about 
5 — 15% of the central r < 2" emission) assuming that 
the X-ray emission from the fifth image has similar in- 
tensity ratios as those measured in the HST ACS image. 
The observed profile may be better fitted by introduc- 
ing the two-component /3-models or increasing the inner 
slope parameter a in the NFW-SSM model; however, we 
will not pursue this process further since the above two 
models already provided acceptable fits to the data. A 
considerably deeper observation is needed to constrain 
further the gas profile in the innermost region. 

From the image analysis, the extent of the X-ray emis- 
sion above the 3-cr background level is found to be 
rx = 91"(= 652 kpc) (Figure El, which is close to an 
overdensity radius of r 50 o — 0.79io.ii Mpc. 

3.3. Cluster Mass Distribution 

Under the assumption of hydrostatic equilibrium, we 
can infer the mass distribution of the lensing cluster from 
the X-ray temperature and surface-brightness profiles. 
Because there is not a significant radial dependence of 
the X-ray temperature in the observations, we assume 
isothermality of the gas in the mass estimation. Figure 



shows the cylindrical cluster mass projected within a ra- 
dius r derived from each surface mass distribution profile. 
We find that both models yield consistent mass profiles 
within the measurement errors. 

It is interesting to compare the mass derived from 
X-ray analysis with th at from gravitational lensing. 
Willi am's fc Sahal l|2004ft derived the cylin drical mass to 
be M (< 100 kpc) = (5± 1) x 1O 13 M , and lSharon eTaO 
(2005) estimated M(< 110 kpc) = 6 x 10 13 M Q . From 
Figure |21 it is clear that both mass estimations are in 
excellent agreement with the mass profile from the X- 
ray observation. For instance, we find the cylindrical 
masses within 100 kpc for the /3-model and NFW-SSM 
model are 5.2± 2 ° x 10 13 M© and 5.0±J;f x 10 13 M Q , re- 
spectively. The discrepancy between X-ray and lensing 
masses has been rep orted in many lensing clusters (e.g., 
lHattori et al.l I1999J) . and is often ascribed to the pro- 
jection of the extra matter along line-of-sight, the elon- 
gation of the lensing cluster along the line-of-sight di- 
rection, or a departure from equilibrium. The excellent 
agreement between X-ray and lensing masses implies that 
none of these effects is significant for SDSS J1004+4112. 

For the NFW-SSM m odel and A c = 187r 2 fi 427 
(Nak amura fc Sutol 11997(1 . the virial mass and the con- 
centration parameter are constrained to be Af V i r = 
6.0±ll x 10 14 M o and c vir = 6.1±i;|. The value of the 
concentration parameter is slightly larger than the theo- 
retically expected m edian value for this virial mass and 
redshift, c vir ~ 4.0 {Bullock et all 12001). but is within 
2-cr scatter among different clusters. 

4. QUASAR IMAGES 

4.1. Anomalous Flux Ratios 

The observed source counts for the four images, A-D, 
are 1237, 1580, 1312, and 763, in the 0.5-7 keV band, 
respectively. Light curves were produced with time res- 
olutions of 2048 or 4096 sec. We did not find any clear 
sign of time variability within the statistical uncertain- 
ties. Note that each component was extracted using a 
circular region with a radius of 2" in the analysis pre- 
sented in this section. The sum of the background and 
contamination from the cluster emission are estimated to 
be only ~ 1% of the source counts for A-C, and ~ 3% for 
D, and are negligibly small compared to the statistical 
errors of the source spectra. 

We measure the energy fluxes of lensed images in the 
Chandra data by fitting each spectrum with a power-law 
plus Gaussian line model (see the next subsection for 
details). The results are summarized in Tabled The X- 
ray flux rati os differ sign ificantly from those measured in 
the optical l|Inada et al.ll2005l) . Since for fold lenses such 
as SDSS J1004+4112, the two images near the critical 
curve are usually brighter than the other two images and 
since the optical flux ratios were r eproduced with sim- 
ple mass models l)Qguri et al.ll2004T . we believe that the 
X-ray, rather than the optical, ratios are "anomalous" . 
The optical/X-ray flux ratios of images C and D are al- 
most the same, therefore the most natural interpretation 
is that the image A is demagnified by a factor of ~ 3 in 
X-rays. This interpretation is consistent with the obser- 
vation that image A appears to be a sadd le-point image 
ijQguri et al.l2004HWilliams fc Sahal2004ft which is more 
likely to be demagnified by perturbation s than the other 
images {Schechter fc Wambsganssll2002|) . While the op- 
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tica l data compared in Table n were taken on 2004 April 

28 iflnada et al] 120051) . the demagnification of image A 
in X-rays is further supported by optical observations on 
2004 Nove mber 27 and 2 005 April 26 with the MAGNUM 
telescope ijYoshii 1 12002). By linearly interpolating these 
two i?-band imaging observations, we infer optical flux 
ratios on 2005 January 1 to be B/A=0.69, C/A=0.40, 
and D/A=0.29, which are not very different from those 
listed in Table □ 

In Table ^ we a l so show the values of a ox to check 
the consistency of the X-ray and optical fluxes. The val- 
ues are somewhat higher than the a — —1.33 that is 
expected from the UV luminosity of SDSS J1004+4112, 
/(2500 A) = 10 29 [ergs -1 Hz -1 ] (assuming a magn ifica- 
tion factor of \i = 50 for image A: lQguri et all 2004) and 
the best-fit a ox -UV lumi nosity relation llSteffen et all 
I2006|) . From Table 5 of iSteffen et~aTl (|2006f) . ~ x is 
-1.408 ± 0.165 and -1.322 ± 0.192 for log £(2500 A) = 

29 - 30 and log £(2500 A) = 28 - 29, respectively. Thus 
the derived a ox is within the l-a (2-cr) scatter for A (B- 
D). 

4.2. X-ray Spectrum 

The spectrum of the lensed quasar component is dis- 
played in Figure 0] To enhance the signal-to-noise ra- 
tio, the total spectrum (sum of the images A, B, C, and 
D) is plotted. We fit the spectrum with a power-law 
model. Again, the Galactic hydrogen column density is 
fixed to TVh = 1-13 x 10 20 cm -2 since we did not find sig- 
nificant intrinsic absorption. This model is rejected at 
the 90% confidence level (x 2 /dof = 144.0/115). We plot 
the residual of the data against the best-fit model in Fig- 
ure which shows ~ 2 — 3d excesses in adjacent bins 
around 2.3 keV in our frame. Then we mask these bins 
and those to either side (9 bins from 2.26-2.66 keV) and 
fit a power-law model. The chi-square becomes signifi- 
cantly smaller (% 2 /dof = 99.5/106). The above results 
suggest the presence of an emission line component at 
~ 2.3 keV. We then fit the spectrum with a power- law 
plus Gaussian line profile model (see Figure^). In com- 
parison to just a power-law model, the fit was improved 
at the > 99.99% confidence level according to the -F-test 
(A% 2 = 35.9 for three additional parameters). The fit- 
ting yields the power-law index T = 1.90 ± 0.04, the cen- 
troid energy of the line E — 6.32^014 keV, and the Gaus- 
sian width of the line a = 355l 12 ;g eV. The line centroid 
is consistent with the neutral iron Ka line at 6.4 keV in 
the quasar rest frame and the derived a suggests that 
the line is intrinsically broad. From the result, we can 
compute the equivalent width EW = 768^ 23 >7 eV in the 
quasar rest frame. The luminosity is 4.8 x 10 45 ergs -1 in 
the 2-10 keV band. In the case of ^ = 50 for the image A, 
though the lens models t hat are consistent with the data 
allow a wide range of \x l)Qguri et alJl2004l) . the intrinsic 
quasar luminosity is estimated to be ~ 2.4 x 10 43 ergs -1 
(2-10 keV). This is within the lu minosity range of Seyfert 
galaxies fe.g.. lGreen et alJl992[) . The reduced chi-square 
of the best fit model is y?j&oi =108.1/112. 

There is a well-known negative correlation between the 
strength of t he iron line and th e X-ray luminosity in 
quasars (e.g., iNandra et al.lll997l) . The derived equiv- 
alent width is larger by a factor of ~ 3 t han t hat ex- 
pected from the relation of INandra et al J l)1997f) and a 



factor of > 4 than that from recent XMM-Newton ob - 
servations l|Page et al.l2004tl.Iimenez-Bail6n et al.l2 005). 
Combined with the result in M4.ll we speculate that mi- 
crolensing demagnifies the X-ray continuum of image A, 
while keeping (or amplifying) the Fe Ka line. More- 
over, microlensing of emission lines in t he optical band 
has b een detected for this lens system ijRichards et alJ 
2004), providing additional evidence for the microlens- 
ing hypothesis. The relative enhancement of Fe Ka line 
has ofte n been observed in the X-ray emission of lensed 
quasars <lOshima et al J 120011: iChartas et all 120021 12004 
iDai et al.ll2003D . and interpreted in terms of microlens- 
ing. Indeed, the relative enhancement of Fe Ka line can 
be re produced if the emission is separated in two regions 
(e.g.. lPopovic et al.ll200l| . 

The origin of the strong Fe Ka line is further explored 
by analyzing the spectrum of each image, rather than 
the sum of all images. To measure the line intensities, 
we fit them with the power-law plus Gaussian model. 
The centroid energy and the width of the line as well 
as the power-law index are allowed to vary in the fit. 
The best-fit parameters and the 90% errors are listed 
in Table We found that the spectral parameters for 
the four images agree with each other at the 95% con- 
fidence level, although the parameters, particularly the 
line width, show large uncertainties. The fitting results 
are also statistically consistent with those derived for the 
total spectrum. We then performed spectral fitting with 
E, a and V fixed to the best-fit values for the total spec- 
trum (Figure 03) and calculated the equivalent widths 
(Table [TJ. Although the errors are very large, the Fe 
Ka line is detected most significantly for image A, sup- 
porting the interpretation above. We also find that the 
hardness ratios of the images are marginally consistent 
w ith eac h other (see Table 

iGreenl l)2006l) recently proposed that spectral differ- 
ences between lensed quasar image components are due 
to small line-of-sight differences through quasar disk 
wind outflows. In particular, the author suggested that 
SDSS J1004+4112 components B-D suffer from absorp- 
tion in a disk-wind outflow and t hat the blue enhance- 
ments in SDSS J1004+4112 A ijRichards et all l2004|) 
could be due to the alleviation of absorption along that 
sightline. Indeed, the A-B difference sp ectrum from 
Keck /LRIS spectroscopy of C IV region ijQguri et alJ 
|2004[) resembles that of the troughs seen in a broad ab- 
sorption line (BAL) quasar. However, the observed line 
profiles and transitions involved (e.g., He II) mean that 
these features more plausibly result from excess emis- 
sion in component A rather than absorption in compo- 
nent B. In addition, quasars with intrinsic UV absorp- 
tion (e.g., BAL troughs) are generally X-ray faint due to 
a large intrinsic X-ray absorption (Nw ~ (0.1 — 10) x 
10 23 cm -2 :e.g. iBrandt et al] l20"00l iGreen et al J 12001 
Gallagher et aljr2002l 12006). however, intrinsic absorp- 
tion was not significantly detected in the X-ray spec- 
tra of the four images; the 90% upper limit on the col- 
umn density is obtained to be N-r < 5 x 10 21 cm -2 , 5 x 
10 21 cm -2 , 4 x 10 21 cm -2 , and 8 x 10 21 cm -2 for A-D, 
respectively. Here we assumed neutral absorbing mate- 
rial at z = 1.734. In the case that the absorber is highly 
ionize d and has high velocity flow, as proposed bv lGreeri 
( 2006) , X-ray spectra with considerably better statistics 
are required to test the model. Moreover, if the differ- 
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ential absorption is responsible for the diff e rent X -ray 
spectra of A and B, we expect ( and iGreenl (|200fif) pre- 
dicts) that the image A should be less absorbed in the 
X-ray. Instead, however, component A has the smallest 
X-ray to optical flux ratio and the X-ray flux ratio is 
A/B = 0.79 ± 0.03 (the 1-er error; see Table QJ. Thus, 
given multiple independent indications, we conclude that 
the microlensing phenomenon is the most natural inter- 
pretation for the present data. 

5. SUMMARY 

We have presented results from Chandra observations 
of SDSS J1004+4112, a large-separation gravitational 
lens system created by a cluster of galaxies. We have 
detected X-ray emission from the lensing cluster as well 
as four lensed quasar images. 

From the cluster X-ray emission, we have constrained 
the bolometric luminosity and the temperature to be 
L x = 4.7 x 10 44 ergs" 1 and kT = 6.4±f;| keV, consis- 



tent with the luminosity-temperature relation of distant 
clusters. We have reconstructed the mass profile of the 
lensing cluster assuming isothermality and hydrostatic 
equilibrium, and found that the mass within 100 kpc ex- 
cellently agrees with that expected from strong lensing. 

X-ray emission from the lensed quasar images displays 
two anomalies: the presence of a strong Fe Ka line and 
significant differences in the flux ratios from those found 
in the optical band. Both of these features are suggestive 
of microlensing. The idea is supported by the fact that 
image A, which appears most anomalous, is a saddle- 
point image. 
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TABLE 1 

X-ray and Optical Properties of SDSS J1004+4112 



Name 




hardness 1 " 


Fo P t c 


Fopt/F-x 




EW [eV] d 


A 
B 
C 
D 


0.95 ±0.03 
1.21 ±0.03 
1.01 ±0.03 
0.56 ±0.02 


0.30 ±0.02 
0.26 ±0.02 
0.27 ±0.02 
0.34 ±0.03 


1.00 
0.73 
0.35 
0.21 


1.00 
0.58 
0.33 
0.35 


-1.27 ±0.01 
-1.16 ± 0.01 
-1.12 ± 0.01 
-1.10 ±0.01 


iici+374 
lioi -554 
crqq+284 
J»J_446 

842_ 389 
t-,-9+369 
00z -268 



a Absorption-corrected, 0.5—7 keV X-ray flux in units of 10 13 ergs 'cm 2 and the 
1-cr error. 

b Hardness ratios defined by the count rate ratios in 2 — 7 keV and 0.5 — 2 keV 
bands. Errors indicate 68% confidence limits. 

c Flu x normalized by the flux of image A in the HST F814W image llnada et alj 
2005). Measurement errors are negligibly small. 

d Equivalent width of the Fe Ka line in the quasar rest frame. 



TABLE 2 

Results of spectral fitting for quasar components 



Name T a E [keV] b a [keV] c EW [eV] d x 2 / dof A X 2 



A 


1 Q2+ - 10 


6 15+ - 26 
D - lo -0.20 


36+ - 17 

U.OD_q 2g 


12521™ 


31.0/20 


18.5 


B 


1.97 ±0.08 


6 45+ 015 


< 0.63 




22.3/20 


9.6 


C 


1.97 ±0.09 


6 25+ ' 35 


0.37(< 0.76) 


1077+J 22 


14.7/20 


14.9 


D 


1.77 ±0.11 


6 12+ 031 
6 S2 +0 15 


0.26(< 0.62) 




30.8/20 


7.3 


A+B+C+D 


1.90 ±0.04 


0.36l u ;i 3 7 


7681^ 


108.1/112 


35.9 



a Power-law index. 

b Centroid energy in the quasar rest frame. 
c Line width in the quasar rest frame. 

d Equivalent width of the Fe Ka line in the quasar rest frame. 
c Chi-square and the degree of freedom of the fit. 

'improvement of the chi-square in comparison to the power-law model. 
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Fig. 1 — Adaptively-smoothed ACIS-S3 image of SDSS J1004+4112 in the 0.5-7 keV band. Both multiple images of SDSS J1004+4112, 
A— D, and the extended emission from the lensing cluster (whose X-ray peak is marked with the cross) are clearly seen. The images A and 
B, which have a minimum angular separation of 3'.'8, are also resolved in the raw image. The point source at (10:04:34.294, +41:12:20.22) 
is a quasar (see text). 




Fig. 2. — Results of the X-ray surface brightness profile fitting with (a) the /3-model and (b) the NFW-SSM model. In each panel, the 
crosses show the observed surface brightness in the 0.5-5 keV band and the solid line shows the best-fit model. The background is shown 
with the horizontal dotted line. The vertical dashed line shows the extent of the diffuse X-ray emission, rx (see text). 
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i 1 1 — i — i — i — i i | 1 1 1 — i — i — i — i t | 




j i i i i_u i i I i Li i i i I i L i I 



0.01 0.1 1 

r [h 70 _1 Mpc] 

Fig. 3. — Enclosed mass of the lensing cluster, Mx, for the /3-model (black) and the NFW-SSM model (red). The dotted lines indicate 
the 90 % error ranges. No te th at Mx is a cylindr ical cluster mass projected within a radius r. The masses derived from gravitational 
lensing I Williams & Saha 200l |Sh aron etail 2005) are also shown for comparison. The meaning of the vertical dashed line is the same as 
Figure |21 
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Fig. 5. — X-ray spectra of images A - D fit with a power-law (dashed) plus Gaussian (dot-dashed) model, which is shown by the solid 
line. 



